Non-gaussian Noise in Quantum Spin Glasses and Interacting Two-level Systems 
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We study a general model for non-gaussian I// noise based on an infinite range quantum fsing 
spin system in the paramagnetic state, or equivalently, interacting two-level classical fluctuators. 
We identify a dilatation interaction term in the dynamics which survives the thermodynamic limit 
and circumvents the central limit theorem to produce non-gaussian noise even when the equilibrium 
distribution is that of non-interacting spins. The resulting second spectrum ('noise of the noise') 
itself has a universal 1// form which we analyze within a dynamical mean field approximation. 
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This classical dynamics is equivalent to the imaginary 
time dynamics of a quantum spin Hamiltonian 



G, 



PACS numbers: 75.50.Lk, 73.50.Td, 61.20.L, 61.20.Lcc 



Noise measurements have long been an important 
probe of the still poorly understood dynamics of spin 
glasses pj , structural glasses and twodevel systems ^,|) . 
Interest in 1// noise |4|,|| has been further enhanced with 
recent progress in the development of nano-scale elec- 
tronic devices such as the single electron transistor ||] 
and precision charge pumps |7j which are extremely sen- 
sitive to offset charge fluctuations B. Such noise also 
plays a significant role in determining the coherence time 
of existing ion-trap realizations of quantum bits || and is 
likely to be even more important in other proposed solid 
state realizations of qubits such as quantum dots {Tot , 
Cooper pair boxes jll], and Josephson junctions 12 1. A 
constant loss tangent in a capacitor leads to a 1// noise 
spectrum |l3) and loss tangents even as a small as 10~ 5 
would have serious consequences in charge-based quan- 
tum bits. Ubiquitous non-equilibrium noise may also be 
implicated in the coherence time of electrons in wires and 
films observed via weak localization and via mesoscopic 
persistent currents | Q . 

We focus here on non-gaussian fluctuations which are 
usually connected with critical phenomena [jig] or sys- 
tems having only a few degrees of freedom j||7 We de- 
velop a simple model with a large number of degrees of 
freedom and identify interaction terms that circumvent 
the central limit theorem and produce non-gaussian fluc- 
tuations in the thermodynamic limit even in the para- 
magnetic state far above any critical point. We analyze 
the non-gaussian fluctuations using the second spectrum 
('noise of the noise') |l(||l|,p"7| . In addition to being of in- 
terest for the study of two-level systems, the dynamical 
model we present is also relevant to models of quantum 
spin glass behavior JlS| . 

Consider TV interacting two-level fluctuators (TLFs) 
labeled by Ising spin variables Si — ±1 whose flipping 
rates obey 
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with the Ising variables S replaced by Pauli spin matri- 
ces a z . For a spin glass, hHi would typically represent a 
uniform external magnetic field, i.e. Hi = 1. For a struc- 
tural glass the hHi represent random local site energy 
differences for the TLFs. jTij is the spin-spin interac- 
tion energy and A = Aj + -^X^jAAj + gG^Sj is 
the barrier height which fluctuates around its bare value 
Aj due to coupling between the TLFs. We take Hi, Aj, 



Gji and IVy 



Tji to be random and uniformly dis- 



tributed on [—1,1], as in the infinite-range Sherrington- 
Kirkpatrick model |Q . 

In steady state, detailed balance gives us the single- 
spin probability ratio 
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Note that this is completely independent of Di which 
controls various features of the dynamics but has no effect 
on the equilibrium statistical mechanics. 

If and only if T is symmetric, we can write the full equi- 
librium probability distribution in terms of a Boltzmann 
factor derived from a Hamiltonian 
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where Z is the partition function. Notice that the matrix 
G need not be symmetric. The approach to equilibrium 
might depend on this, but the equilibrium does not. 

For large enough 7 (low enough temperature) the sys- 
tem undergoes a glass transition into a frozen state. 
Near this transition the dynamics may also become non- 
gaussian, but we restrict our attention here to the high 
temperature paramagnetic regime where 7 is small. Our 
central result is that even if 7 = so that the equilibrium 
distribution is that of non-interacting spins, the second 
spectrum indicates non-gaussian dynamics for any A 7^ 0. 
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The term Yljjti ^jSj ^ s a 'dilatation' term. It causes 
all the barriers to collectively and simultaneously increase 
or decrease in a correlated manner. In contrast to this, 
interactions through the gGij and 7^ terms cause bar- 
rier fluctuations which are statistically independent on 
different sites (at least for small g and 7 so that the sys- 
tem is above the critical point). As a consequence of the 
central limit theorem, one expects that only the dilata- 
tion term survives to produce non-gaussian fluctuations 
in the thermodynamic limit of this infinite-range model. 

A possible physical realization of such a dilatation term 
is the elastic interaction between TLFs. In addition to 
the dipole term one could have an expansion of the lat- 
tice coupled to the impurity positions [fl9|| . The motion 
of a defect then changes the hydrostatic pressure by an 
amount proportional to ^ throughout the volume V . 
Since the barriers would tend to go up uniformly with 
increasing hydrostatic pressure, a long-range dilatation 
term with an amplitude ~ A -1 would be present [fl9|| . 
Such a term thus could be important in small grains. We 
have chosen for simplicity to scale the interaction with 
A - ' in order to obtain a sensible thermodynamic limit 
in the present model. 

To ensure that the power spectrum of the model decays 
as 1/f, we take the bare barriers to have a flat distribu- 
tion on a fixed interval [Q A l G [A min , A max ] . We carry 
out simulations using Eq. (|]J). To average, each "mea- 
surement" was repeated 1600 times, where A», Aj, Gy 
and Tij are assigned new values each time. Using St as 
a discretizing time interval, the noise signal is defined as 
V(t St) = N-y^f^Siit St). Henceforth, we set St = 1. 
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FIG. 1. Power spectrum as a function of frequency for a 
non-interacting (A = g = h = ~f = 0) and two interacting 
(A = g = h = 7 = 0.25, 0.5) cases. Other parameters are: 
A" = 14, A min = 1, M = 2 22 and N = 100. (W{v)) is 
found to be independent of N . Thin dotted line is a plot of 
the analytical expression for the non-interacting case, Eq. (Q) . 
Typical barrier fluctuations for the (X = g — h = ^y = 0.5) 
case are ~10% and ~100% of A max and A min , respectively, 
indicating weak or moderate interactions. 



In this letter, t, ti and v, are dimensionless integer 
times and frequencies, respectively. The normalization 
A -1 / 2 ensures that (V 2 {t)) = 1, independent of A. 

To study the time-time correlation of the noise, 
we use the power spectrum defined as (W{v)) = 

(|M- 1 / 2 £fl 1 e 2 ™' t / M F(t)| 2 ), where M is the length 
of the time sequence. We see in Fig. [l] that the power 
spectrum is insensitive to weak interactions. 

To study non-gaussian properties of the noise, we use 
the second spectrum introduced by Restle et al. |T(||l],|l7| . 
Given a time sequence of length M, we split it into Mi 
shorter time sequences of length Mi each, and define the 
second spectrum as 
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The quantity Wt 2 (y) is a single shot "power spectrum" 
of V(t) in the time window t%. Physically, the second 
spectrum is a spectrum of the power spectrum's wander- 
ing within a chosen frequency band, [z/ mm , i/ max ]. It is 
a four-point correlation function which probes the noise 
of the noise [|l6|,|l],[l7j . For the non-interacting (gaussian) 
case the second spectrum becomes flat Q and takes the 
value E^=7m.n (W{v)) 2 . We see in Fig. | that for A = 
the second spectrum decreases for increasing A, merging 
into its non-interacting (gaussian) counter part for large 
N. However, for A = 0.5 the second spectrum is inde- 
pendent of A for large A. Thus, the dilatation term in 
Eq. ([|) is the only term that survives the central limit 
theorem. 
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FIG. 2. Second spectrum as a function of frequency for 
A = and 0.5 with N = 100, 200, 400, 800, 1200. Other pa- 
rameters are: g = fe = 7 = 0.5, Mi = M 2 = 2 10 , A max = 15, 
A min = 0.5, z/ min = 125 and ^ max = 175. Data for the A = 
case is shifted down one decade for clarity. Dots (0) show the 
second spectrum for the non-interacting (gaussian) case. 
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All the other terms produce only gaussian fluctuations 
in the thermodynamic limit. 

We now calculate the power spectrum and second spec- 
trum analytically. First, we set g — h = 7 = since, 
as shown above, these terms do not contribute to 11011- 
gaussian properties of the model and the power spec- 
trum is insensitive to weak interactions. Second, we 
set Ai = (A max - A min )i/N + A min . Third, we re- 
place Aj € [— 1,1] by Ai = 1. Mathematically, it is al- 
ways possible to make a variable ('gauge') transforma- 
tion Sj = sign(Aj)Sj such that all Aj are positive in the 
new basis. Simulations have shown that the results for 
these two cases are qualitatively the same, with the sec- 
ond spectrum of the Aj = 1 case somewhat larger as 
expected. 

We now use the Suzuki- Trotter decomposition [|TJ and 
recast the quantum model, Eq. (||) , into a classical model 
with the partition function 



Z = 



exp 
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Here, V(t) is an auxiliary field used to decouple the 
spins. Furthermore, e;(i) = jge~ Al ~ A ^ (t) and Ki(t) = 
iln(coth[ei(i)]), where (3 = l/k B T is the inverse 
temperature. In Eq. (^), we have replaced S[V(t) — 
N- 1/2 E^S,(t)} by S[V(t)-N-^j:f =1 S 3 (t)} and, 
thereby, introduced an error of order 0(N~ X / 2 ) which 
vanishes in the N — > 00 limit. 

Using Eq. (^|), the power spectrum for the non- 
interacting case may in the large N and large M limit 
be shown to take the well known form 



(w>;A = 0)) = 



atan (lFF7M m) ~ atan ( e ™/M m) 



(7) 



(A max - A min )(irv/M) 

From Eq. (m), we see that for intermediate frequencies 
(W(v;X — 0)) ~ 1/l>, independent of the temperature. 
Henceforth, we set /3/M — 1 so that the analytical and 
simulation results are directly comparable. 

We plot Eq. (0) in Fig. [l] using the same parameters 
as in the simulations. There is nice agreement between 
the simulation and analytical results. At the highest fre- 
quencies, v > 10 6 , results from the simulations are higher 
compared to the analytical result. This is expected since 
in the simulations the time sequences are discrete, while 
in the analytical calculations a continuum description is 
used to integrate out the sums. It is precisely at these 
highest frequencies that the difference between the two 
methods becomes important. 

We now apply the dynamical mean-field approximation 
(DMFA) jl8| by replacing the formal averaging over V in 



Eq. (JsJ) with the requirement that all (multi-point) cor- 
relation functions of V(t) are self-consistently identical 
to the corresponding correlations functions of V(t). For 
example, (V(t)V(t')) v = (V(t)V(t')) s . Here, (> g and 
Qy are averages over the spin configurations and the 
auxiliary field, respectively. As for other infinite-range 
spin models p8j , this approximation becomes exact in 
the N — > 00 limit. 

For a classical ID Ising model (S(t)S(t')S(t")S(t'")) = 
(S(t)S(f)) (S(t")S(t"')}, given that the overlap between 
the time intervals [t,t'] and \t" ,t'"\ is zero. Thus, the 
term {W t2 {v) W t '(i^)) s in Eq. ^ may be replaced 

by (WtM) s (Wt> 2 ("')) s 



0(M 2 1 ), where the error 



1/2 

0(M 2 1 ) vanishes in the M2 — > 00 limit. Using Eq. 
and the DMFA we obtain 
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We split V into fast and slow components, V = Vf + V s , 
and then assume that V s is constant over the time inter- 
val [t', t"] whereas J2t=t' e~ XVf ^ is replaced by its second 
cumulant average. The second cumulant depends only 
on the first spectrum, which from Fig. (|l|) we know is 
universal and independent of the all couplings including 
A. In Eq. (0) we may therefore make the approxima- 



tion E i=t ' 



-\V(t) 



\t' - t"| e - AV '=( t 2)+ 2 2-« where £ is an 



unimportant constant which will not contribute to the 
second spectrum computed to order A 2 below. Applying 
the approximation to Eq. (B) we find 
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within the frequency band 1 -C v A A . Here, 

v e [0, ^]. Eq. (§) shows that the shape of (W t2 {v; X)) s 
is not modified by the interactions between the fast and 
the slow fluctuators. However, its amplitude shifts up or 
down depending on the configuration of the slow TLFs 
during the short time interval ti- In other words, the 
wandering of the power spectrum is controlled by the sig- 
nal of the slow fluctuators. Because the barrier distribu- 
tion is flat, there is no net effect on the first spectrum, but 
the mean square wandering produces the second spec- 
trum. 

Expanding Eq. (|^) to second order in A and ensemble 
averaging over the fluctuations of V s (t,2) we see that the 
non-gaussian contribution to the second spectrum is, up 
to a simple normalization factor, determined by the slow 
part of the first spectrum: 
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From Eq. @ and 
the universality of the first spectrum shown in Fig. (|l|) 
we see that, apart from the prefactor, the shape of the 
non-gaussian part of second spectrum is independent of 
all interaction parameters as well as the size and posi- 
tion of the frequency window used to probe it, provided 
we restrict the frequency band to be inside 1 <C v <C 
^e~ A ~ A . In Fig. || we plot simulation results for 
four different sets of parameters. We see that the sim- 
ulation results collapse onto the line predicted by the 
analytical result. The small deviations may be explained 
by the fact that A is not infinitesimal, by practical lim- 
itations on the available dynamical range which cause 
the inequality below Eq. (||) to be imperfectly satisfied, 
and by the unavoidable poor sampling statistics for the 
slowest fluctuators. 

In the present model, the second spectrum has been 
found to decay close to 1//; that is, the 1// noise it- 
self has 1// noise. This is the same result seen experi- 
mentally in semiconductors [^2],f|]|] , discontinuous metal 
films |p2fl , spin glasses Q and carbon resistors [jl7|. Our 
model is an infinite-range interaction model in which 
we have chosen the simplest possible dynamics that cir- 
cumvents the central limit theorem to yield non-gaussian 
noise even in the thermodynamic limit of the weakly cor- 
related paramagnetic phase. To the extent that the Ising 
variables we use could represent entire droplets of spins, 
our model could also potentially be relevant to the 
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FIG. 3. Second spectrum minus it gaussian counterpart as 
a function of frequency for g = h = 7 = 0, A max = 17, 



A mm = 1) ^ = M% = 2 L* &nA jy = 10Q ^ A W (2) {u 2 ) ) 

is found to be independent of N for N > 100 when 
g = h = 7 = 0. Solid line is a plot of the analytical ex- 
pression Eq. (|lo|). 



strongly correlated regime of spin glasses Q . 
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